Single quasiparticle excitation dynamics on a superconducting island 
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We investigate single quasiparticle excitation dynamics on a small superconducting aluminum 
island connected to normal metallic leads by tunnel junctions. We find the island to be free of 
excitations within the measurement resolution allowing us to determine Cooper pair breaking rate to 
be less than 3 kHz. By tuning the Coulomb energy of the island to have an odd number of electrons, 
one of them remains unpaired. We detect it by measuring its relaxation rate via tunneling. By 
injecting electrons with a periodic gate voltage, we probe electron-phonon interaction and relaxation 
down to a single quasiparticle excitation pair, with a measured recombination rate of 8 kHz. Our 
experiment yields a strong test of BCS-theory in aluminum as the results are consistent with it 
without free parameters. 



The quasiparticle excitations describing the micro- 
scopic degrees of freedom in superconductors freeze out 
at low temperatures, provided no energy exceeding the 
superconducting gap A is available. Early experiments 
on these excitations were performed typically close to the 
critical temperature with large structures so that Ns, the 
number of quasiparticle excitations, was 

hightlhi- 

Later 

on, as the fabrication techniques progressed, it became 
possible to bring Ns close to unity to reveal theparity 
effect of electrons on a superconducting island !10l4l4[. 
In recent years, the tunneling and relaxation dynamics 
of quasiparticles, which we address in this letter, have 
become a topical subject because of their influence on 
practically all superconducting circuits in the low tem- 
perature limit |15l422l |. 

We study the quasiparticle excitations in a small alu- 
minum island shown in Fig. [1] (a). The island is con- 
nected via a thin insulating aluminum oxide layer to two 
normal metallic copper leads to form a single-electron 
transistor (SET) allowing quasiparticle tunneling. We 
bias the SET by a voltage VJ, between the source and 
drain and polarize the island by an offset charge Ug ex- 
pressed in units of e with a gate voltage Vg ~ erig/Cg. 
Here Cg is the gate-island capacitance. The current / 
through the SET is governed by sequential tunneling of 
single quasiparticles and it exhibits Coulomb diamonds 
which overlap each other because of the superconducting 
energy gap 2^ 23] , observed for our structure as a region 
bounded by the red sawtooths in Fig.[T](b). 

In the sub-gap regime, eVb < 2 A, of Fig.[T](b) the cur- 
rent should be suppressed if there are no quasiparticle 
excitations present. Nonetheless, we observe a finite cur- 
rent which has a period twice as long in ng as compared to 
the high bias region, a unique feature of a superconduct- 
ing island due to Cooper pairing of electrons. Its origin 



is a single electron unable to pair in the condensate and 
hence remaining as an excitation. This parity effect has 
been observed in the past in similar structures lol- 12 
but typically with two-electron Andreev tunneling be- 
ing the main transport process. We focus on devices 
where Andreev current is suppressed since a high charg- 
ing energy, Ec > A, makes tunneling of two electrons 
energetically unfavourable compared to that of a single 
quasiparticle [23J26j. In this case, the transport is dom- 
inated by single-electron processes allowing simple and 
direct probing of the quasiparticle excitations without 
the interfering multi-electron tunneling. 

For a quantitative descripition of the transport char- 
acteristics, we performed a numerical simulation of the 
device operation shown in Fig. [T] (c). To describe simul- 
taneously the charging of the island with electrons and 
the excitations involved in superconducting state, we as- 
sign probability P{N, Ns) for having N excess electrons 
and Ns quasiparticle excitations on the island. The time 
evolution of P{N, Ns) is described by a master equation 



PiN,Ns)= J2 

N',N' 



where P{N, Ns) stands for the time derivate oi P{N, Ns) 
and rjv'-s.jv,Af^-*-Ars for the transition rate from state 
P{N',N'g) to P{N,Ns). These rates are set by electron 
tunneling between the island and the leads. Cooper pair 
breaking and recombination of quasiparticles. 

Tunneling rates are calculated by the standard first 
order perturbation theory so that electrons tunneling into 
the superconducting island to a state with energy E > 
A will increase the quasiparticle number and electron 
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FIG. 1: (color online), (a) Scanning electron micrograph of 
the sample studied. It is biased with voltage T4 and a gate 
offset voltage Vg is applied to a gate electrode. The latter is 
not shown in the micrograph. The operation is governed by 
few quasiparticle excitations which can be either particle or 
hole like. A particle like excitation can relax by tunneling out 
(filled circle) and a hole by electron tunneling in and filling 
the state, (b) Measured source-drain current J as a func- 
tion of bias and gate voltages, (c) Calculated current based 
on sequential single-electron tunneling model, (d) Measured 
current at rig = is shown as black dots. Black line is calcu- 
lated assuming a quasiparticle generation rate of 3 kHz, and 
red one assuming a vanishing generation rate. 



number by one, 



Tn- 



e^Rr 



dEns{E){l-fs{E))fN{E + SE), 

(2) 

where Rt is the tunneling resistance, ns{E) is the BCS 
density of states, fs/N{E) the occupation probability of 
state E in superconductor / normal metal and 6E the 
energy gain from charging and biasing [2^, [23]. If the 
incoming electron tunnels to the lovifer branch E < —A, 
it removes an excitation by filling a hole. Similarly, an 
outgoing electron from the upper branch removes an ex- 
citation while if it tunnels out from the lower branch, an 
excitation is generated. 

In the tunneling rates the exact quasiparticle number 
Ns is accounted for by a nonequilibrium distribution fs ■ 
In general this has a complicated form as a function of 
energy. However due to the fact that quasiparticles are 
injected close to the gap, the resulting tunneling and re- 
combination are not sensitive to the functional form of 
fs- Our system is also symmetric with respect to the two 
branches. Therefore we neglect branch imbalance and 
paramerize the quasiparticle number by an effectively in- 
creased temperature Ts and a Fermi distribution in case 



of fs- This gives the relation 

Ns = V2^D{EF)VVAkBTse-^^''^^' , (3) 

where D{Ef) = 1.45 x lO^'^J-^m^^ is the density of 
states in the normal state [i^l and V the volume of the 
island. For the normal metallic leads we use Fermi distri- 
bution with Tjv = 60 mK, equal to the base temperature 
of the cryostat. For a detailed description of all transition 
rates, see the supplementary information [2^ . 

For the steady-state represented by Fig.[T](c), we solve 
Eq. ([T|) with P{N, Ns) = and calculate the current as 
an average of the tunneling rates weighted by the prob- 
abilities P{N,Ns)- The parameter values of sample A: 
Ec = 240 /ieV, A = 210 /leV and tunneling resistances 
Rti — 220 kfi and Rt2 = 150 kft for the two junctions 
were used in the simulations. They were determined from 
measurements in the high bias regime (e|T4| > 2 A) and 
hence their values are independent of the sub-gap fea- 
tures. 

The simulation of Fig. [T] (c) reproduces the behaviour 
observed in the experiments. The relaxation rate of a 
single quasiparticle excitation via tunneling is expected 



to be Fa 



{2e^RTD{EF)Vy 



190 kHz, where we have used the measured dimensions 
for V = 1.06 /im x 145 nm x 25nm. From the fit in the 
sub-gap regime we obtain Fqp = 150 kHz, consistent with 
the prediction. The value of Fqp affects only the value of 
current on the light blue plateau of Fig.[T](c), not the ac- 
tual form or size of the terrace. As our simulation based 
on sequential tunneling reproduces all the features in the 
sub-gap regime, we conclude that the two-electron peri- 
odicity originates from single-electron tunneling and the 
operation is essentially free of multi-electron tunneling 
processes. At odd integer values of rig, the characteristic 
feature of Andreev tunneling would be a linear-in-Vb cur- 
rent at low bias voltages and a subsequent drop 11, 2^, 
which is absent in our data. The leakage current in the 
sub-gap region does not vanish even in the zero temper- 
ature limit but remains essentially the same as presented 
in Fig. [TJ Therefore, all quasiparticle excitations can- 
not be suppressed at finite bias voltages by lowering the 
temperature, if Ug is close to an odd integer. 

At even integer values of Ug we have ideally no cur- 
rent flow as all electrons are paired. If Cooper pair 
breaking would take place in the island, we would ob- 
tain two quasiparticle excitations in the superconductor. 
One of these excitations can then relax by tunneling to 
the leads followed by tunneling of a new excitation to 
neutralize the offset charge, similarly as at the odd Ug 
case described above. This cycle would continue until 
the two quasiparticles recombine to a Cooper pair. As 
the recombination rate, discussed below, is slower than 
the rates in the cycle, we have several electrons tunneling 
through the device for each broken pair, hence amplifying 
the signal. With this model we obtain an upper bound 
^ N^Nfl^2 = 3 kHz for the pair breaking rate under our 
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experimental conditions based on simulations shown in 
Fig. [l] (d). This generation rate corresponds to energy 
absorption at 2rjv^Ar.o^2^ ==0.1 aW power on the su- 
perconducting island. 

Under constant biasing conditions, there is at most 
one quasiparticle present in the subgap regime at low 
temperatures. The non-tunneling relaxation on the is- 
land is then not possible, since recombination would call 
for two excitations. Therefore the static case can be de- 
scribed by pure tunneling without other relaxation pro- 
cesses. To study recombination of two quasiparticles into 
a Cooper pair, we injected intentionally more quasipar- 
ticles to the island. The injection was done by a periodic 
drive of the gate voltage. By changing Ug, we change the 
potential of the superconducting island and either pull 
quasiparticle excitations into the island when the poten- 
tial is lowered or create hole type excitations as potential 
is raised and quasiparticles tunnel out. The number of 
injected quasiparticles and the number of quasiparticles 
on the island can then be determined from the resulting 
current curves with the help of simulations. In the ex- 
periment we approach two different limits which we will 
discuss in the following: When the pumping frequency 
is high, Ns is large. Such situation can be described by 
a thermal models i.e., by an increased time-independent 
effective temperature of the superconducting island. In 
the opposite limit of low frequency, Ns is small. Then 
the thermal model fails and we have to account for the 
exact time dependent number of quasiparticles. 

In Fig. [2] (a) we show the measured current for three 
different values of bias voltage Vf,. The gate drive is sinu- 
soidal around Ug = 1/2 with amplitude Ag expressed in 
units of e/Cg and frequency / = 1 MHz corresponding 
to fast pumping. Without accumulation of quasiparticles 
to the island, the current would show quantized plateaus 
with spacing e/, similar to the SINIS turnstile [23|. How- 
ever, as the island has a surplus of quasiparticles, i.e. it 
is heated up, the current is substantially higher. We 
use now the thermal model where the heat injection to 
the island by electron tunneling is balanced by electron- 
phonon interaction. The heat flux into the phonon bath 
is given by 

■ de {nie, Ts) - n(6, Tp)) /_°°^ dE 



24C(5)fc^ 



xnsiE)ns{E- 



ifsiE) - fsiE + e)) . 



(4) 

where E is the material constant for electron-phonon cou- 
pling, (^(z) the Riemann zeta function and n{t,T) = 
(exp(e/(fcsr)) — the Bose-Einstein distribution of 
the phonons at temperature Tp. See supplementary in- 
formation for derivation. The same result is obtained by 
kinetic Boltzmann equation calculations [sol , 31 1. The 



simulations based on the thermal model are shown as 
black lines in Fig. [2] (a). The electron-phonon coupling 
constant E — 1.8 • 10^ WK~^m~^, used in simula- 
tions, was measured in the normal state, where Qop = 
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FIG. 2: (color online), (a)-(b) Measured current of Sample 
A against gate voltage amplitude Ag ai f — 4 and 1000 kHz 
at bias voltage values Vb = 120, 200, 280 /iV shown as blue, 
red and green dots respectively. Black lines show simulations 
assuming an elevated temperature on the superconducting is- 
land, (c) Similar measurement for Sample B at / = 1 MHz. 
(d) The electron-phonon heat flux in the superconducting 
state normahzed by that in the normal state extracted from 
the measurements (black circles). Temperature is expressed 
with respect to the critical temperature Tc = A/1.76A;b. The 
theoretical result of Eq. Q is shown by the black line. Solid 
and dotted red lines show the recombination and scattering 
part of Eq. Q correspondinly. The open grey symbol is from 
Sample B. 



T.V{T^ - T|) and Tn is the electron temperature 

We expect the thermal model to be a good approach 
if Ns ^ 1. With the high frequency and large ampli- 
tude drive in Fig. [5] (a), we have Ns ^ 10 quasiparticles 
present for Ag ^ 1, suggesting that the thermal model is 
adequate for these data. If the frequency is lowered to 
/ = 4 kHz, shown in Fig. [21(b), the thermal model fails as 
Ns approaches unity. As a further proof of the overheat- 
ing, we repeated the high frequency measurement using 
sample B with measured parameter values Ec = 620 /xeV, 
A = 270 ^leY, Rn = 1800 kfl, Rt2 = 960 kfl and 
V — 800 nm x 60 nm x 15 nm. The result is shown in 
Fig[2](c). Again, the simulations (black lines) are able to 
reproduce all non-trivial features of the measured curves. 
As a summary of the thermal model fits, we repeated the 
measurement of Fig. [5] (a) at different frequencies and 
determined by numerical simulations the temperature of 
the superconducting island and the heat injected into it 
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FIG. 3: (color online). Red circles show the measured cur- 
rent for / = 4 - 1000 kHz at Vt = 280 fiY. Black lines 
are simulations based on Eq. ([l]). Recombination rates are 
taken to match the heat flux in the regime where the thermal 
model applies (Eq. (01). Dotted blue lines are calculated with 
vanishing electron-phonon relaxation rate and solid gray lines 
with the thermal model. 



based on the measured current. The results are shown in 
Fig. [2] (d) as black circles. The results match well with 
the expected electron-phonon coupling of the supercon- 
ductor, Eq. (HI), presented as the solid black line. 

For studying the deviation from the thermal model, 
we measured the characteristics at four frequencies / = 
4, 10, 100 and 1000 kHz at Vb = 280 /iV, shown in Fig. H 
The thermal model is presented now as solid gray lines. 
A more adequate description of the data at low frequen- 
cies is obtained by simulations based on Eq. ([T]), where 
we keep track on the number of excitations during the 
operation and take into account the recombination rates. 
The heat flux of Eq. (|4|) consists of recombination and 
scattering which can, respectively, be expressed as 



Qs 



ttVT, 
305)11" 



{kBTsA^ 

-A/keTs 



|(fcBT5)2A3)e-2A/'=BT.^ 



(5) 

Their contributions are presented in Fig. [2](d). Whereas 
the scattering does not change the number of quasiparti- 
cles Ns, recombination leads to transitions Ns — Ns — 2. 
We account for this process by including the recom- 
bination rate rjv-s.Af,A's-^JVs-2 — Qrec(-^s)/2A, where 
the relation between the effective temperature Ts in 
Eq. ([5]) and the exact quasiparticle number Ns is given 
by Eq. The solid black lines of Fig. [3] are calcu- 

lated with the same value of E as obtained in the normal 
state. Blue dotted lines show similar simulations where 



electron-phonon relaxation is disregarded. 

With the simulations based on instantaneous quasi- 
particle number we can reproduce the experimental fea- 
tures precisely with no free parameters in the calcula- 
tion. At the lowest frequency, / = 4 kHz, we have only 
one quasiparticle present for most of the time. Hence 
the curves are not sensitive to the recombination. As 
frequency is increased, the simulations without electron- 
phonon relaxation deviate from the experimental data. 
For / = 10 kHz we probe the recombination rate of a 
single qp pair only, Tn^n^2^o = 8 kHz. We checked 
this by artificially changing the recombination rates for 
Ns > 2, without any significant difference in the curves. 
At higher frequencies, the recombination for Ns > 2 be- 
comes significant as well. The results of the two mod- 
els approach each other and the thermal model becomes 
valid. 

In summary, a small superconducting island at low 
temperatures has allowed us to study the dynamics of 
single electronic excitations and their relaxation. Under 
quiescent conditions we found a vanishing Cooper pair 
breaking rate within the measurement resolution: based 
on the measurement noise we obtained an upper limit of 
3 kHz for this rate. On the other hand, by periodically 
pumping electrons, we controllably increased the number 
of quasiparticles and were able to measure the recombi- 
nation rates both in the large quasiparticle number limit 
and for a single quasiparticle pair: rjv_i.Ar,2^.o — 8 kHz. 
The recombination rates are in quantitative agreement 
with the relaxation measured at higher temperatures and 
in the normal state. 
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MASTER EQUATION 

We model the superconducting island of the main article by writing down a master equation by assigning a probability 
P{N, Ns) for having N excess electrons and Ns quasiparticles on the island. The time derivate of P{N, Ns) is then 



P(7V, Ns) 



N'. N' 



r 



N' 



.N,N',^NsPiN',N's), 



(1) 



where Tn'^n,n' -^Ns is a transition rate from state {N',Ng) to {N,Ns). For N' — N,Ng = Ns we insert a rate 



N^N,Ns^Ns 



N'N' ^ N^N'^N,Ns- 



^N'g=iNsj which is a sum of all rates out from the state (N,Ns). Next we 
consider the rates arising from electron tunneling between normal metallic leads and the island and the recombination 
and pair breaking in the superconductor. 



ELECTRON TUNNELING 

Within first order, the tunneling rates are given by Fermi golden rule similarly as in Refs. 1^. Now as we consider also 
explicitly Ns, we need to split the rates so that an electron tunneling into the upper branch of the superconducting 
island to energy E > A, increases the excitation number Ns by one and an electron tunneling into lower branch, 
E < — A, decrease it by one as shown in Fig.[TJ Likewise, an electron tunneling out from the upper branch decreases 
Ng and an electron tunneling out from the lower branch increases it by one. Hence we obtain rates 



- N-l^N, Ns-l^Ns 



N-l^N,Ns + l^Ns 



E 



1 



e2i?T. 



j JA 



E 



1.2 



E- 
E 



dE nsiE)il - fs{E, Ns))fNiE + SEj{N)) 

A 

dE ns{E){l - fs{E, Ns))fN{E + SEj{N)) 
dE ns{E)fsiE, Ns)il - fN{E + 6Ej{N))) 
dE ns{E)fs{E, Ns){l ~ fN{E + SE.iN))), 



-A 



(2) 



where Rtj is the tunneling resistance of junction j, ns{E), the BCS density of states, fs{E,Ns), occupation proba- 
bility for state E on the island if we have Ns quasiparticles. and fN{E) occupation probability in the normal metallic 
lead at energy E. SEj{N) the energy cost from biasing and charging energy We assume /jv to be a Fermi dis- 
tribution with temperature equal to the bath temperature T^r = 60 mK. For the superconductor we need to have fs 
such that the number of excitations equals Ns. For the tunneling rates, as well as for the electron-phonon interaction 
discussed below, the functional form of the distribution is irrelevant as long as the excitations are close to the gap 
edges E « ±A, where we inject them. We also neglect the branch imbalance, i.e. assume no chemical potential shift. 
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Figure 1: (color online). The four different tunneling processes changing A*' and Ns- (a) Electron tunneling out of the island 
and removing an excitation, (b) Electron tunneling into the island creating and excitation, (c) Electron tunneling into the 
island and removing a hole excitation, (d) Electron tunneling out and creating a hole excitation. 

Sfi = 0, because our drive is symmetrical with respect to the two branches. Hence we use a Fermi distribution with 
temperature T5. The relation between Ns and T5 is then 

dE-^===fE^V2^DiEF)V^/Ak^e-^/''^^^, (3) 

where D{Ep) is the normal state density of states at Fermi level and V the volume of the island. The tunneling 
current is obtained then as 

N^N+l, Ns^Ns±l — Ns-^Ns±l ) P ^s) ■ (4) 

N,Ns,± 

Including higher order tunneling to Eqs. ([1]) and (0]) is straightforward. For example Andreev current gives rise to 
terms r]\[±2^N,Ns^Ns- These rates are calculated similarly as in Refs. 0-111 since quasiparticle number Ns is not 
changing. Including Andreev tunneling terms into the calculations in main text, does not change the results as Ec > A. 
For Ec < A, those terms contribute because Andreev tunneling thresholds are exceeded before the single-electron 
tunneling thresholds 0-01 ■ 



ELECTRON-PHONON INTERACTION 

To obtain the recombination rates for Eq. ([1]) and Eq. (4) of the main article we use Hamiltonian 

H = H,+Hp + H,p. (5) 

Here Hg + Hp is the non-perturbated part and Hep the perturbation from electron-phonon interaction. The BCS 
Hamiltonian of the superconductor is 

with Ek = \J e^, + I Afcp. The diagonalizing fermionic operators 7^-0- are connected to the electron operators cwa' by 
cL = Ka^-{ka) + "fca7L' where we use notation 

/ Vfe if (T =t / Uk if cr =t 

The coefficients satisfy |ufcp -f = 1 and /Vf^Vk/uk = Ek — ek leading to 

kP-l-KP-i(l-;|). (8) 



Similarly, we have for the phonons 



Hp^J2f2Wqb%. (9) 
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The coupling of these two systems has the form 



= wy2(4cfc_<,65 + clck+qbl). 

k,q 



as in Ref. p . The operator of heat flux to the phonons is 



Hp = ^[H,Hp] - ^[H,p,Hp] = iiyJ24^\cW-qbq - clck+gbl). 

k,q 



By using the Kubo formula in the interaction picture we obtain 



H„) = (H„) -- 
h 



dt' 



Hp{t), Hcp{t') 



k' a' 



lit) Ck+q„{t) b\{t), cl,^,it') Ck'-qa'{t') bqit') 



(10) 



(11) 



(12) 



Next we use bq{t) = e *"<!*5^ and c\^{t) = vl^-f_(ka){t) + "^kallA^) '^ilh ^kcj{t) = IkaC ^^kt/h_ example, let us 

consider the first non-zero term in Eq. (|12p . We obtain 



(13) 



iu,{t'~t)i(-Ekt+Ek-^t-Eyt'+Ey+^t')/h* * /r t h -y , -y'^ 



ii^,(t'-t) i(Ek-,~Ek){t-t')/h\ |2| |2/r t L t Ltl\ 



where we have paired the operators as fc = fc' + q. For some of the terms, there exists two different pairing possibilities. 
From the second commutator of Eq. (jl2p we find similarly a term 

(14) 



where the dummy summation index has been changed from k' to k. Now we can combine these two terms and do the 
time integration to obtain energy conservation rules as delta functions. By combining all terms similarly, assuming 
Afc = Ae"^, where A is real, and changing the summing to integration we obtain from Eq. ([T^ 



TTV^N{EF)D{q) 







-Q 






EkE^^ 


'Q 






-Bfe Ek- 


-Q 


A^ 




EkEk- 


-Q 



(1 - fk) fk-q (nq + 1) - fk (1 - fk-q) nqj6{Ek " Ek-q + fuVq) + 
(1 - fk) (1 - fk-q) {riq + 1) - fkfk-qnq)S{Ek + Ek-q + ftcj,) + 
' fkfk-q {nq + 1) - (1 - fk) (1 - fk-q) nq]d{-Ek - Ek-q + hujq) + 
Jk (1 - fk-q) (nq + 1) - (1 - fk) fk-qnqjS{-Ek + Ek-q + fiWg)|, 



(15) 



where N{Ef) is the density of states of the normal state and D{q) 

is the volume of the system. Now, / d'^q — I-k dqq^ J^^ d{co8{9)) where 9 is the angle between k and q. Further 



j^yj the density of states of the phonons. V 



more 



2m 



2m 



h^kp 



qcos{9) and ujq — ciq, where q is the speed of sound. For integrating the 



delta functions over cos(6'), we need 



dE, 



k — q 



d{cos{9)) 



E, 



k — q 



ek-q h?kFq 



ns{Ek-q) 



h?kpq ' 



(16) 
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where we have identified the density of states of the superconductor ns{E) = ^ 



E 



Therefore we obtain 



27r 



■N{EF)ci 

A 



dq 



m ( 



EkEk 


-9 


+ 




^ EkEk 


-g 


+ 




^ EkEk 


-g 






EkEk 





) ( (1 - /fe) (n, + 1) - /fe (1 - n,) 
(1 - /fe) (1 - /fe-g) {ng + 1) - /fc/fe-,n, 

/fe/fe-g + !)-(!- fk) (1 - /fe-g) "-q 

) (/fe (1 - fk-g) [ng + !)-(!- fk) fk-gn^) 



Ek-,=Ek+hi^g>0 
Ek-g = -Ek-hujg>0 

Ek-q = -Ek + huJg>0 ^ 

Ek-q=Ek-hu)g>Q > 



(17) 



Now, we change the integration variables to Ek = i/e^ + and e — hug = hciq. Hence, we have dek 
and dq — {hci)~^de. We can also use notation E = Ek and E' — Ek^g. Then, equation p7)) yields 



ns{Ek)dEk 



Qcp 



miy^VN{EF) 



dE 



de e ns 



iE)ns{E'){ 



(l - ( (1 - fm f{E') {n{e) + 1) - f{E) (1 - f{E')) n( 

+ (l + w) ( (1 - fiE)) (1 - f{E')) (n(e) + 1) - f [E] f {E')n{e)) 
+ (l + ^) {f{E)f{E') {n{e) + !)-(!- /(S)) (1 - /(i?')) n{e)) 
f{E) (1 - /(i?')) {n{e) + !)-(!- /(i?)) /(i?')n(e) 



A^ 
EE' 



E'=E+e>0 
E' = -E-e>0 
E' = -E+e>0 
E' = E-e> 



(18) 



Finally we can change integration variable E and the dummy variable E' as follows: for first term: no changes, for 
second term: change of sign for £", for third term: change of sign for E and for the fourth term: change of sign for 
E and E' . Then, we see that both E and E' cover both positive and negative values and we obtain 



rnv'^VNiEp) 
TrhPkpcf 



dE / dee^nsiE)ns{E + e)x 



E{E+c) 

With some algebra, we can simplify Eq. ([T7| to 



(l-TO)) {i^-fiE))fiE + e) inie) + l)-f{E)il-f{E + e))n{e)). 



(19) 



Qcp — 



de (n(e, Tg) - n(e, Tp)) / dE ns{E)ns{E + e) 1 - 



24C(5)fc| 7o 

Here we have identified the electron-phonon coupling constant S 



(/(i?)-/(i? + e)). (20) 



2AC(b)kgmv^N{EF) 



by taking the limit A 



0, 



which yields the normal-state result Qcp — SV(r| — Tp). One arrives to Eq. (|20p also by using kinetic equations for 
the electronic excitations under electron-phonon interaction [7[ . 

For obtaining the recombination and scattering rates, we identify in Eq. (llSp the first and last term to correspond 
to scattering and the two middle ones to recombination and pair breaking. Since the heat rate is set dominantly 
by the quasiparticle occupation near the gap, we can again parametrize the number of quasiparticles by an effective 
temperature Tg and the Fermi distribution fp = {exp (E/k^Ts) + 1)~^ without any loss off generality. For Tp <C 
Ts ^ A/fcs, satisfied in the experiments, we obtain the recombination and scattering heat fluxes as 



l-eciTs) 



= iTs) 



■Kb 

^^DiEp)V 



di 



d^'{E + E') 



i\3 



1 



A2 



EE' 



/eJe 



^(..T.A^ + I(..T.)W)e-/--, 



^DiEp)V 
n 



di 



dCO{E-E'){E-E- 



1 - 



A2 



EE' 



Instead of parametrisation with respect to Ts, we could also parametrize the efi'ective temperature Ts by the quasipar- 
ticle number Ns with Eq. ([3]) and ask for the rate rjv->.Ar.Ars+2-^7Vs entering the master equation for P{N, Ns). This 
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accounts for the transition of Ns + 2 to Ns quasiparticles, whereas the number of excess electrons N is unchanged. 
As the quasiparticles lie close to the gap the energy loss during this process is 2A and therefore the rate is given by 



N^N,Ns+2^Ns 



Q,ec{Ns + 2) 



2A 

In these terms the heat balance then is fully described by the master equation for P{N, Ns). 



(21) 
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